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Local orbits of a pure state of an TV x TV bi-partite quantum system are analyzed. We compute 
their dimensions which depends on the degeneracy of the vector of coefficients arising by the Schmidt 
decomposition. In particular, the generic orbit has 2N 2 — N— 1 dimensions, the set of separable states 
is 4(iV — 1) dimensional, while the manifold of maximally entangled states has TV 2 — 1 dimensions. 
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O ' I. INTRODUCTION 

<n : 

| The existence of entangled states, i.e., roughly speaking, the states of a composite system which exhibit quantum 
correlations among the subsystems, appeared recently to be extremely important in rapidly developing field of quantum 
communication. It is due to non-classical properties of entangled states that various schemes of quantum computing, 
quantum cryptography and quantum teleportation can be thought of being practically realizable. 

A pure state \ip) in the Hilbcrt space H = Ha KHb of a composite quantum system consisting of two subsystems 
A and B with Hilbcrt spaces Ha and Hb is separable, if it can be cast to the product form \tp) = \iPa) <S> \i^b), where 
■ \4>a) and \4>b) are some states of the subsystems. States which are not separable are called entangled. The situation 
is more complicated in the case of a mixed state (a density matrix p) jl| . It is separable if it is expressible as a convex 

sum of product states: p — J2iPiPi^ ® p\ B \ Pi > Oj J2iPi — 1; where p-" 4 ' and p\ B ^ are, in general mixed, states of 
the subsystems. A mixed state is called entangled if it is not separable. In what follows we consider only systems with 
finite-dimensional Hilbert spaces which seem to be more important in proposed applications of quantum information 
,_H \ theory, the infinite-dimensional case needs some refinement of the above definition of separability. 

It is relatively easy to check whether a given pure state is separable or entangled (e.g. by investigating its Schmidt 
coefficients - see below). The situation complicates for mixed states - we do not know how to check unambiguously 
separability of a given mixed states if the dimensionality of the Hilbert spaces of subsystems exceeds 3 j| . 

As a problem complementary to determining the separability properties of a given state one can pose the question 
of the relation between the set of the separable (entangled) states to the set of all states of the composite system. This 
can be understood as the question of a relative measure of the set of entangled states (i.e. "how probable is that a 
given state is entangled?" ) - the problem posed and partially solved in , or about the geometrical and topological 
properties of this set. In this paper we concentrate on the latter problem in the following setting. Since we are 
interested in quantum correlations between two subsystems we should take into consideration only these properties 
. which do not change under various quantum mechanical operations performed locally in each subsystem. Thus two 
states which are interconvertible one to another via local unitary transformations (i.e. purely quantum mechanical 
operations without decoherence) are equivalent from the point of their entanglement properties. This can lead to 
construction of appropriate measures of entanglement characterizing the classes of equivalent states. Our approach 
is in a sense complementary to the task of identifying the set of all invariants with respect to the local unitary 
transformations [p-pd|. 

In this work we pose and solve the question of the dimensionality and topology of manifolds of states equivalent to a 
given one via local unitary transformations. Thus the present paper may be regarded as an extension of ]T^ | (see also 
[p~3| |l5|), in which these questions were discussed for the simplest system of two qubits. In the case of pure states we 
find the explicit results for any N x N composite system by identifying explicitly the topology of the orbits as well as 
in a purely algebraic, algorithmic manner. The second approach which does not depend on the Schmidt decomposition 
(see below) is, in principle, applicable also to mixed states, this is illustrated by considering the generalized Werner 
states flj. 
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II. PURE ENTANGLED STATES 



A. Schmidt decomposition 

Consider a pure state \ip) of a composite Hilbert space TL = TCa'&'Hb of size N 2 . Introducing an orthonormal basis 
{l n )}^Li m each subsystem, we may represent the state as 

N N 
n— 1 m— 1 

The complex matrix of coefficients C of size N needs not to be Hermitian nor normal. Its singular values (i.e. the 
square roots of eigenvalues A& of the positive matrix C^C) determine the Schmidt decomposition |l6|-|lgfl 

N 

W)=Y J V^k\k')®\k"), (2.2) 
fe=i 

where the basis in 7i is transformed by a local unitary transformation W (g> V. Thus \k') = W\k), and \k") — V\k), 
where W and V are the matrices of eigenvectors of C^C and CC^ , respectively. In the generic case of a non- 
degenerate vector A, the Schmidt decomposition is unique up to two unitary diagonal matrices, up to which the 
matrices of eigenvectors W and V are determined. The normalization condition = 1 enforces EitLi = !• 

Thus the vector A = (A l7 Ajv) lives in the (N — 1) dimensional simplex Spj. The Schmidt coefficients A& do not 
depend on the initial basis \n) ® |m), in which the analyzed state is represented. 



B. Pure state entanglement 



The Schmidt coefficient of a pure state are equal to the eigenvalues of the reduced density operator, obtained 
by partial tracing, p A = tr b {\ip) ■ A pure state is called separable, if it can be represented in the product form 
IVO = IVu) ® IV*b)) where \ipA) & Ha an d \^b) G Hb- This occurs if and only if there exists only one non-zero 
Schmidt coefficient, Ai = 1, i.e. the reduced state p A is pure. In the opposite case state is called entangled. 
A pure state is called maximally entangled if all its Schmidt coefficients are equal, Ai = Xk = 1/N. Note that the 
Schmidt coefficients are invariant with respect to any local operations Ul = Ua <8> Ub, and thus they may serve as 
ingredients of any measure of entanglement. 



C. Local orbits 



We are going to study the orbits of a given pure state with respect to the local transformations [/ £ . Two 
states belonging to the same orbit are called interconvertible, since they may be reversibly transformed by local 
transformations one into another fl9[| . Let us order its Schmidt coefficients A = (0 < Ai < A2 < • • • < Ajv). In order 
to describe the character of the degeneracy we rename them into A = (0, • • • , 0, V\, ■ ■ ■ , V\, v%, ■ ■ ■ , vi, . . . , vk, • • • , vk) 
where each value v n occurs m n times and mo is the number of vanishing Schmidt coefficients. Obviously too + 



E 



K 



m r , 



N , and Too might be equal to zero. The main result of our paper is contained in the following 



Proposition. The local orbit generated from has the structure of the following quotient space 



O = 



U{N) x U{N) 
<?(m ,TOi, . . . ,m K )'' 



(2.3) 



where Q(mo,mi, . . . ,rriK) is the subgroup of the direct product U(N) x U(N) consisting of the pairs of unitary 
matrices (U, V) of the form 



Uo 



vo 



U 



V 



(2.4) 



U K 



"K 
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where uq and vq arbitrary matrices from [/(mo), and u\, . . . ,uk denote arbitrary matrices from, respectively, 
U(mi), . . . , U(rriK)- The overall phase factor e 1 ^ accounts fo r th e irrelevant phase of the state \ip), ie. we identify 



states differing by a phase factor. The dimension of the orbit (2.3) reads 



K 



dim(O) = 2N 2 - 2m 2 - Y^ m l ~ 1 ■ ( 2 - 5 ) 



Indeed, let us observe that the action of the tensor product U <S> V £ U(N) ® U(N) on the state (2.1) 



U®V\^)=Y J C mn U\m)®V\n)= £ C mn U km \k) ® V ln \l) = J2(U CV T ) kl \k) ® \l), (2.6) 

m,n m,n,k,l k,l 

reduces to the direct product action on the coefficient matrix C 

U(N) x U(N) 3 (U, 7):Ch (17, V){C) := UCV T (2.7) 

Let now the action of (U, V) reduces C to its diagonal Schmidt form 

UCV^ — diag(0, . . . ,0,1*,... ,v 2 ,... ,v K ,... ,v K ). (2.8) 



Then UCV T = UUCV T V T iff U and V are given by (2.4) Now the formula (|2.3|) follows in an obvious manner, once 



we realize that in fact we should disregard any unimportant overall phase of (|2.l| ) (or in other words we should ide ntify 



the coefficient matrices C and C = e C, ie. work in an appropriate projective space). The dimension formula (2.5) 



follows from a simple calculations involving the dimensionalities of the unitary groups, while the last term equal to 



unity stems from the projectivisation procedure. An alternative, algebraic proof of this result is proved in Section III. 
In fact the orbit has a structure of a Cartesian product: 

= ffi) x U{N) (29) 

C/(m ) x U{mx) x ■ ■ • x U(m K ) U(m ) x U(l) ' v ' ' 

where the first factor represents global orbits in the set of density matrices of size N with the same spectrum [^0|,^lj . 
In the language of fiber bundles such orbits form the base, while the fibers consists of all N x N pure states, which are 
related by partial tracing to a given density matrix of size N. We shall provide a complete proof of this fact elsewhere 
II- 

In the generic case of all coefficients different (and non zero), i.e. K = N the manifold is thus identified as 



_ u(N) um 

Ua ~ p7(l)P X [/(!)' 



with the dimension 

dim(C 3 ) = 2N 2 - N - 1 . (2.11) 

The set of all orbits enumerated above produces the complex projective space <CP n2 ~ 1 - the (2iV 2 — 2) d imensional 
manifold of pure states of the N x N system. However, the the set constructed of the generic orbits ( ]2.11 ) generated 



by each point of the interior of the Weyl chamber, is of full measure in the space of pure states. In this way we 
demonstrated a foliation of CP N _1 . This foliation is singular, since there exist also (measure zero) leaves of various 
dimensions and topology, as listed in Table 1 for N = 2, 3 and 4. 

D. Special cases: separable and maximally entangled states 



For separable states there exists only one non zero coefficient, Ai = 1, so mo = N — 1. Thus (2.3) gives 

U(N) U(N) 

u(i) x u(n-i) x up) x uysr - 1) 



O scP - TT ^ U i N i T\ x ^ = CP— x CP—, (2.12) 
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with the dimension dim(O scp ) = 4(N — 1). The maximally entangled states are characterized by Ai = An = 
hence mi = N and mo = 0. Therefore 

_UW_SU(N) 

with the dimension dim((9 max ) = iV 2 — 1. Note that this space is not isomorphic with SU(N) because U(N) is not a 
direct product of [7(1) and SU(N) f||. Since SU(N) x U{1) = U(N) x Z N , where Z N is the discrete permutation 
group of N elements, the orbit of the maximally entangled states can be written as O max = SU(N)/Zpf. This 
structure follows also from the fact that the entire orbit may be written as max = (U <8> 1)1^), where \ty) is an 
arbitrary maximally entangled state, and U is an arbitrary unitary matrix determined up to an overall phase 63]. 



E. Special cases: N = 2, 3 and 4 



The set of all possible Schmidt vectors A form the N — 1 dimensional simplex Sn- Its corners represent N mutually 
orthogonal separable states, while its center denotes the maximally entangled state = Q^fcLi \kk))/y/N. Any 
permutation of the Schmidt coefficients may be obtained by a local transformation of the pure state. Therefore it 
is sufficient to consider the orbits generated by Schmidt vectors belonging to a certain asymmetric part <Sjv of the 
simplex, so called Weyl chamber. Any ordering of the Schmidt coefficients corresponds to choosing one chamber out 
of Nl, in which the simplex Sn can be decomposed. 

The Schmidt simplex and exemplary Weyl chamber for N — 2,3 and 4 are presented in Fig.l. (Note that the 
simplex of diagonal density matrices of size N, obtained from N x N pure states by partial tracing, has the same 
geometry). The numbers by each part of the boundary of Sn denote the dimensions of the local orbits, which are 
listed in Table 1. In the simplest case N = 2 the simplex reduces to the interval [0, 1], while its asymmetric part S2 
equals to [0,1/2]. The edge generates the four dimensional orbit of separable states, CP 1 x CP 1 , and the point 
1/2 leads to the 3-D manifold of maximally entangled states max = SU(2)/Z2 ~ £0(3) ss RP 3 . This structure was 
pointed out by Vollbrecht and Wern er |2^] , and the above singular foliation of CP 3 was discussed in |l^-[]l|. In the 
case of any point inside the simplex fl2.12|) gives the following topology of the generic 2x2 local orbit 



°> = uW x W) = 3 x MP ' (2 - 14) 

in agreement with recent results of Mosseri and Dandoloff JIB] . 



III. ALGEBRAIC DETERMINATION OF ORBIT DIMENSION 



A. General case: N x N mixed states 



The reasoning presented in the previous section hinges on the Schmidt decomposition of the density matrix for a 
pure state. As such it cannot be extended to mixed states. For this reason we present an alternative method introduced 
in |12[ ], which can be, in principle, applied also in the latter situation. It is based on purely algebraic reasoning, and, 
as such, gives only local information, i.e. only about the dimensions of the manifolds of interconvertible states and 
not about their topology. 

Although the group of local unitary transformations is C — U(N) ® U(N), it is obvious that since its elements act 
on an arbitrary density matrix p £ C N <g> C N by conjugations, p 1— * UpU\ we can take in fact C = SU(N) ® SU(N) 
instead. Let R 2 ^ 2 - 1 ) 3 s ^ U(s) £ SU(N) <g> SU(N) be some parameterization of the group C such that U(0) = I 
(i.e. s = (si, S2, • ■ ■ , S2JV 2 -2) are the coordinates in SU(N) with the origin at the unit matrix). The tangent space 
to the local orbit through p (i.e. to the space of the states interconvertible with p) at this point is spanned by the 
vectors: 

p u := Af/( s)pt /t (s) | s=o . (3.1) 

The dimension of the tangent space, hence of the manifold itself, equals the number of linearly independent vectors 
Pk- 
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From the unitarity of U(s) it follows: 



Pk 



dU 



,P 



s=0 



= [h,p\ = pl- 



(3.2) 



The number of independent pk equals the rank of the 2(N 2 — 1) X 2(N 2 — 1) Gram matrix (the unimportant factor 
of 1/4 is introduced for further convenience) 



Gran • — ^^PmPr, 



which, upon using (|3.2|), can be cast into 



Choosing the standard parameterization of SU (N) in the vicinity of the identity we obtain 

mr 

9s k j s _„ 



k ■- 



ie k (g)I, k = 1, . . . ,N 2 - 1 
il®e k , k = N 2 ,... ,2N 2 -2, 



where ejt = — et are generators of the Lie algebra su(N). They obey the commutation relations 

[ej,ek] = Cjkiei, 

where Cjki denote the structure constants and we use the summation convention. We normalize e k to fulfill 

Tr BjBk = —2Sjk- 

An arbitrary hermitian matrix p acting in £g> C N can be decomposed with the help of su(N) generators 

p ■= jr^I + iak{ek ® I) + ik(I <g> ei) + C mn (e m ® e„). 



From (gj), (p.5|). (pO|), (p77p, and (gj) the Gram matrix |D|) is calculated as 

G 



4 B 



(3.3) 

(3.4) 

(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 



where the (TV 2 — 1) x (TV 2 — 1) matrices A, B, and D read 



Amn = C m jkC n lk(N CLjdl + 2Cj r C/ r )/2, 
Bmn = C-mjkC-nlrCklC jr ■) 

(Nbjbi + 2C rj C r i)/2. 



(3.10) 



B. Special case: N x N pure states 

Using the abo ve o utlined procedure we can recover the results for pure states obtained in Section I. For a pure state 
p = \tp) {ip\ Eq. ( pTll ) reduces to 

G 

ran — 

W^W- (3.H) 



We choose the following explicit form of the generators expressed in the standard basis {|1) , |2) , • • • , |iV)} of C 

/ fc \ 



^yy ^|fc+l)(fe+l|-^|0(/|J , fc=l,2,...,JV-l, (3.12) 
= (m\ + \m) (n|), 1 < m < n < N, (3.13) 



fc(fc 

e$ n = \n) (m| - |m) (n| , 1 < m < n < N. (3.14) 
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We reorder the non-diagonal generators e>mn and e„ by changing two indices {mn} into a single one k according to 
k = N — 1 + (to — 1)A — to(to + l)/2 + n in the case of and k = N - 1 + A(A - l)/2 — to(to + l)/2 + n in the 
case of emli, so that {e^}, k = 1, 2, . . . , TV 2 — 1 is the desired complete set of generators. 

It proves to be more convenient to use not Ik themselves, but the following linear combinations of them: 



L h :=i(e fc ®7 + Z®e fc )/2, 1 < k < A 2 - 1, 
£fc :=i(e fc ® J- I<g>e fc )/2, N 2 < k < 2A 2 — 2, 



(3.15) 
(3.16) 



what amounts to a mere change of basis in the Lie algebra and, obviously, does not influence the rank of G. 

After rather straightforward but lengthy calculation we find G in the form ( |3.9| ) with 5 = and bloc-diagonal 
matrices A and D 



A = 



A™ 

A( 2 ) 



D = 



£>« 
£>( 2 ) 



(3.17) 



The blocks A^ 2 ) and Z)( 2 ) are diagonal (A 2 — A) x (A 2 — A) matrices with the diagonal entries 

1 < k < (A 2 - A)/2, 



A {2) - 



N 



3=1 



4*=(V^-V^) 2 E A i ' (A 2 — A)/2 < k < N 2 



A, 



vi =1 



(2) 
fcfc 



(VA 



%f E A ? > 1 < fc < (A 2 - A)/2, 



vi =1 

A? 



(2) 
fefc 



(0^+ Va7) 2 [ ] > (A 2 — A)/2 < k < N 2 



A, 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



In each of the above formulas (to, n) is the unique pair of numbers such that < to < n < A and fulfilling 
(to - 1)A - to(to + l)/2 + n = k for 1 < k < (A 2 - A)/2 or (to - 1)A - m(m + l)/2 + n = k - (A 2 - A)/2 for 
(A 2 - A)/2 < k < A 2 - A. Moreover, we find that of two (A - 1) x (A - 1) matrices and £)W the latter equals 
zero, while the former reads 



A {1) = 



(Efcli Afc _ wA m+i ) ( X)feLi Afc ) - (Efcli Afe - m\ m+1 ) (Ek=i Afe - iA„+i) 



\/ to(to + l)n(n + 1) 

1} (Efe=l A fc + n2X n+l) (EftLl A *) - (Efe=l A fc ~ «A„+l) 2 



nm 5 



m < n, 



A {1) = 



n(n + 1) 



(3.22) 



(3.23) 



A 



In this way we found that the entire matrix G has at least A — 1 vanishing eigenvalues (due to £)W = 0), A 2 
doubly degenerate eigenvalues (A; ± A.,) 2 (the eigenvalues of A^ and D^) and the A — 1 eigenvalues of A^\ 

Although, at first sight, A^ looks quite complicated, it is relatively easy to calculate the traces of its powers 
Tr(AW) fe , k = 1, 2, . . . , N — 1 and, consequently, its characteristic polynomial 



N 



k=l 



P (A) := det{A^ - A) = ^ (-l) fe+1 k Pk X N - k . 
Here pi = n, p 2 = r 2 , and p fe = r fc fe^Li • 



(3.24) 



fc-2 



where Tfc are the coefficients of 



N 



Q (A) := J] ( A - A = E T " AAr ^'' 



(3.25) 



k=l 
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i.e. the elementary symmetric polynomials in Ai, A2, . . . , Xn of the order k. Observe that due to the normalization 



Sfe=i Afc = 1 we can substitute Tk for pk in (3.24) and, consequently. 



N 



P(\)=Y^ k Pk X N - k = \Q> (A) - NQ (A) , (3.26) 

k=l 

where Q' (A) := dQ (X)/dX. It follows immediately that the multiplicity of the root A = in P equals the multiplicity 
of A = in Q (i.e. the number of Schmidt coefficients equal to 0). Indeed, if Q(X) = X k Qi(X) and Qi(0) 7^ then 
Q'(X) = fcA fc - 1 Qi(A) + A fe Q' 1 (A) and P(A) = X k [{k - N)Qi(X) + XQ[(X)} = A fc Pi(A), where P r (0) = (fc- AT)Qi(O) ^ 
since k < N. 

Now we are ready to calculate the rank of G. There are 

1. N — 1 vanishing eigenvalues of D^ 1 ' , 

2. too vanishing eigenvalues of A^\ 

3. for each m K -de generate Schmidt coefficient to„(to„ — 1) vanishing eigenvalues of of the form ( |3.19 ) and of 
the form j3T2C|) of P/ 2 >, 

4. 2too(toq — 1) vanishing eigenvalues of and of the forms ( |3 . 1 8| |3~21 ). 



hence the co-rank (the number of zero eigenvalues ofG) equals (N — 1) + mo + J2n=i(. m n ~ m «) + 2(to,q — to ) = 2m 2 ,- 
J2n=i m n"lj where we used m + Yl n =i TO ™ = N. Consequently, taking in account that G is an 2(N 2 — 1) x 2(N 2 — 1) 



matrix, its rank equal to the dimension of the orbit is given by (2.5). 

As mentioned at the beginning of the section, the above analysis can be, in principle, extended to mixed states. To 
show this let's consider (admittedly rather trivial) example of the generalized Werner states 

p=i^J + t#>M, (3-27) 

where the pur e state l^) is characterized by the Schmidt numbers (0 < Ai < A2 < • • • < Xn). It is obvious that the 
Pk of Eq. ( p.2[ ) are, up to the scaling factor a the same as for the pure state \tp). Consequently, the dimension of the 
orbit through p is determined by the Schmidt coefficients of \ip) exactly in the same way as previously. 



IV. COEFFICIENTS OF THE CHARACTERISTIC POLYNOMIALS AS ENTANGLEMENT MEASURES 



There exist several non equivalent ways to quantify quantum entanglement [)25|-|27|. Following Vedral and Plenio 
p8f we assume that any entanglement measure 

i) equals to zero for any separable state, 

ii) is invariant with respect to local unitary operations, 

iii) cannot increase under operations involving local measurements and classical communication. 

For pure states, p = these requirements are fulfilled by the Shannon entropy of the Schmidt vector, (in 

other words von Neumann entropy of the partially reduced density matrix), Ei(\ip)) = — 53fe=i Afc m ^fc> simply called 
entropy of entanglement, as well as the generalized Renyi entropies, E a (\ijj)) = in(J2k=i ^fc)/(l — a ) p9|^C 



Consi der n ow the coefficients Tk of the characteristic polynomial (3.24) of the nontrivial block A^ 1 ' of the Gram 



matrix ( 3.11 ) for a pure state of a N x N bipartite system. As derived above they are given by the elementary 



symmetric polynomials in Ai, A2, • • ■ , Ajv of the order k 

N 

ri = 5^A fc = l, 

k=l 
N N 



t 2 = y,J2 A * A '' 

k=l l=k+l 
N N N 

k=l l=k+l m=l+l 
N 

T N =n x k . (4.1) 

fe=i 
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Due to the definition of the Gram matrix the coefficients Tfe, k = 2, . . . ,N are invariant with respect to local unitary 
transformations and are equal to zero if and only if the state is separable. 

As shown recently by Nielsen any pure state \tp) may be transformed locally into a given state \<j>), if and 
only if the corresponding vectors of the Schmidt coefficients satisfy the following majorization relation -< A^. Any 
entanglement measure cannot increase under such an operation. This condition is fulfilled by the coefficients Tfe, 
since the elementary symmetric polynomials are known to be Schu.r-conca.ve functions J32| , for which A -< jl induces 
T (^) > T (^)- Thus the quantities (4.1) posses the property of entanglement monotones, and their set consisting of 
N — 1 independent elements, {t2, . . . tn}, provides the complete characterization of the pure states entanglement p9[ . 
Beside the simplest case of N = 2, (for which all measures of the entanglement generate the same order in the set of 
pure states p0| ) , the coefficients are not functions of the Renyi entropies and induce different orders in the set of 
pure entangled states. 

It might be interesting to analyze how the traces of the Gram matrix, t% := tr(G fe ), change during non- unitary 
local transformations. Our numerical experiments performed for mixed states of 2 x 2 system suggest that all traces 
tk, k = 1, . . . , 6 do not increase under local bistochastic transformations , p t-y p' = YliPiU^ <8> Uf pllf^ ® Uf\ with 
'YliiVi = 1- The question whether this property holds also for systems of higher dimensions remains open. 
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FIG. 1. Simplex of Schmidt coefficients Sn for pure states of N x N system with N = 2, 3, and 4; (the same picture may 
also represent the set of the spectra of density matrices of size iV obtained from pure states by partial tracing). Right hand 
side shows an asymmetric part of Sn - the Weyl chamber Sn, while the numbers denote the dimensionality of local orbits 
generated by each point. 
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Table 1. Topological structure of local orbits of the N x N pure states generated by one Weyl 
chamber of the simplex of the Schmidt coefficients, D s is the dimension of the subspace, while D 
represents the dimension of the orbit, ( o ) denotes separable states, while ( * ) denotes maximally 
entangled states. 
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